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We update our calculations of thermal-photon production in nuclear collisions at the Relativis- 
tic Heavy-Ion Collider (RHIC). Specifically, we address the recent experimental observation of an 
elliptic flow of direct photons comparable in magnitude to that of pions, which is at variance with 
expectations based on quark-gluon plasma (QGP) dominated photon radiation. Our thermal emis- 
sion rate is based on previous work, i.e., resummed leading-order QGP emission and in-medium 
hadronic rates in the confined phase. These rates are nearly degenerate at temperatures close to the 
expected QCD-phase change. The rates are convoluted over an improved elliptic-fireball expansion 
with transverse- and elliptic-flow fields quantitatively constrained by empirical light- and strange- 
hadron spectra. The resulting direct-photon spectra in central Au-Au collisions are characterized by 
hadron-dominated emission up to transverse momenta of ~ 2-3 GeV. The associated large elliptic 
flow in the hadronic phase mitigates the discrepancy with the measured photon--^ compared to 
scenarios with QGP-dominated emission. 



I. INTRODUCTION 



Dileptons and photons are the only particles which es- 
cape the interior of the fireball created in ultrarelativistic 
heavy-ion collisions (URHICs) unaffected. Their produc- 
tion, on the other hand, is rooted in the strongly inter- 
acting medium and thus illuminates the latter's prop- 
erties, see Refs. [IHH for reviews. Recent highlights of 
electromagnetic observables in URHICs include "mea- 
surements" of the in-medium vector-spectral function via 
dilepton invariant-mass spectra at the SPS and of 
the medium temperature via direct photons at RHIC Q ■ 

It is well known that the thermal component in the ob- 
served electromagnetic spectra results from a convolution 
of the temperature- and density-dependent emission rate 
over the entire space-time history of the expanding fire- 
ball in URHICs. From the interplay of decreasing tem- 
perature and increasing three- volume in the course of the 
fireball expansion, it can be fairly well established that 
photon and dilepton emission in the low-energy regime, 
5o ^ 1 GeV 1 , are dominated by the hadronic phase. At 
energies beyond ~ 1 GeV, the situation is less clear, 
since the competition between hadronic and QGP sources 
will be sensitive to additional ingredients [101 ] . e.g., the 
relative strength of hadronic and QGP emission rates, 



the phase-transition temperature (which formally demar- 
cates the space-time dependence of the two sources) and 
the transverse flow (inducing a blue shift to higher q t 
which is more pronounced in the later hadronic phase) 2 . 

In this context, recent measurements of the elliptic flow 
of direct photons (i.e., after the subtraction of long-lived 
hadron decays, mostly 7r°, rj — > 77) in semicentral Au-Au 
collisions [llj have revealed remarkable results. It has 
been found that the pertinent flow coefficient, v^qt), 
is as large as that of charged pions up to momenta of 
q t ~ 3 GeV (albeit the photon data carry somewhat 
larger error bars). This result is difficult to reconcile with 
a dominant QGP emission source. Primordial photons 
from binary N-N collisions, whose emission is expected 
to be isotropic, will further reduce the total direct-photon 
t>2- Current model calculations (l2H14| using a hydrody- 
namically expanding medium with QGP and hadronic 
radiation, as well as primordial photons, underprcdict 
the experimentally measured V2 by a factor of ~ 5 (~ 3 
when accounting for the maximal systematic error in the 
measurement). The question thus arises what could be 
missing in these calculations. 

In the present paper we re-examine several aspects 
related to thermal photon emission from the hadronic 
medium. Since the V2 in the hadronic phase of URHICs 
is large, an augmented hadronic component is a natural 



1 The energy variable, qg, encompasses both mass (M) and 2 Note that the blue-shift distortion does not apply to dilepton 
transverse-momentum (qt) dependencies, e.g., qo = (M 2 +5 2 ) 1 ' 2 invariant-mass spectra if a finite detector acceptance can be fully 

for q z = 0. corrected for. 
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candidate to improve the description of thermal-photon 
emission at RHIC and thus reduce the discrepancy with 
the V2 data. First, we note that typical hydrodynamical 
evolutions with first-order phase transitions tend to un- 
derestimate the radial (and possibly elliptic) flow built 
up in the fireball at the end of the QGP phase. This was 
borne out of a recent phenomenological analysis of light- 
and strange-hadron spectra using blast-wave parameteri- 
zations of the respective sources at thermal and chemical 
and freezeout (with T ch ~ 100 MeV and T fo ~ 180 MeV, 
respectively) [15j . In particular, the observed universal- 
ity in kinetic-energy scaling of the vi of light and strange 
hadrons, together with an earlier decoupling of multi- 
strange hadrons (as inferred from their p t spectra), sug- 
gests that most of the hadronic V2 is indeed of partonic 
origin [jtjj . Second, for the thermal emission rates from 
hadronic matter we use our previous results of Ref. @, 
which, in particular, include contributions from baryons 
which arc known to be important from dilepton calcu- 
lations [lfj, even at RHIC [T^]. Finally, we take into 
account chemical off-equilibrium effects in the hadronic 
phase, i.e., effective chemical potentials for pions, kaons, 
etc., which can further augment the hadronic compo- 
nent in thermal-photon spectra (e.g., typical processes 
like 7rp — ¥ 7T7 are enhanced by a pion fugacity to the 
third power). 

Our paper is organized as follows. In Scc.[II]wc briefly 
recall our input for the thermal-photon rates as taken 
from Ref. as well as for the primordial contribution, 
which we check against pp data. In Sec. lIIII we update our 
description of the thermal fireball evolution by construct- 
ing (a time evolution of) flow fields which is consistent 
with the empirical extraction at chemical and thermal 
freezeout. In Sec. IIVI we evolve the rates over the fire- 
ball evolution and discuss the resulting photon-g t spec- 
tra, in particular the composition of the thermal yields. 
In Sec. |V] we present our results for the direct-photon 
V2 in comparison to the recent PHENIX data [ll[ and in 
light of other model results. Finally, Sec. IVII contains our 
conclusions. 



II. PHOTON SOURCES 



The photon spectrum resulting from a heavy-ion (or 
pp) reaction is usually referred to as inclusive photons. 
The subtraction of long-lived final-state decays leads to 
the notion of direct photons, which are the ones of inter- 
est in the present context. For direct photons, we further 
distinguish the radiation of thermal photons fSec. Ill Af . 
characterized by an equilibrium-emission rate to be inte- 
grated over the space-time evolution of the medium, and 
a non-thermal component emanating from primordial in- 
teractions prior to thermalization fScc. Ill Bl) . 



A. Thermal Emission 

In the present work we adopt the thermal emission 
rates of photons as developed and compiled in Ref. Q. 

For QGP radiation we use the numerical parameteri- 
zation of the complete leading-order in a s rate as given 
in Ref. [l8|]. The main input for the QGP rate is the 
strong coupling "constant" for which we take an expres- 
sion with temperature-dependent one-loop running at the 
scale - 2ttT, a s (T) = 67r/271n(T[GeV]/0.022). This 
amounts to values of around 0.3 in the relevant tempera- 
ture regime, T = 1-2T C , which turns out to be consistent 
with recent estimates from the (pcrturbative ) C oulomb 
term in in- medium heavy-quark free energies [22( . 

The basis of the thermal emission rate in hadronic mat- 
ter forms the electromagnetic correlation function com- 
puted in Refs. [23[ using hadronic many-body the- 
ory with effective Lagrangians. It has been successfully 
used In in the interpretation of dilepton data at the 
SPS [J, M and has been carried to the photon point 
in Ref. [9|. It includes a rather extensive set of meson 
and baryon resonances in the interaction of the isovec- 
tor current with a thermal heat bath of hadrons. It 
has been augmented by additional meson-exchange re- 
actions in a meson gas which become important at pho- 
ton momenta q > 1 GeV, e.g., n, lo, and a\ exchange in 
7T + /3— >7r + 7as well as strangeness-bearing reactions 
(e.g., 7r + K* —> K + 7) An important element in 
constraining the hadronic vertices to empirical informa- 
tion, such as hadronic and radiative decay branchings, 
is the (gauge-invariant) introduction of vertex form fac- 
tors. The latter lead to a substantial reduction of the 
hadronic emission rate with increasing photon momen- 
tum, which is essential for quantitative descriptions of 
hadronic emission rates at the momenta of experimental 
interest (qt < 3.5 GeV for thermal radiation). Without 
vertex form factors, hadronic photon rates should not be 
considered reliable for momenta q > 1 GeV. We do not 
include here additional tttt Bremsstrahlung contributions 
(mr ititj) as evaluated in Refs. [2^.[25|. This source is 
important at low momenta, qt < 0.5 GeV, put plays no 
role in the region of interest of the present investigation 
(q t > 1 GeV). 

The resulting total hadronic and QGP emission rates 
were found to be remarkably close to each other for tem- 
peratures around the putative transition temperature, 
T c ~ 180 MeV. This is appealing from a conceptual 
point of view in terms of a possibly continuous matching 
of the bottom-up and top-down extrapolated hadronic 
and QGP rates, respectively. It is also welcome from a 
practical point of view, since it much reduces the uncer- 
tainties associated with identifying the medium in the 
fireball as hadronic or partonic, which, in the case of a 
cross-over, may not even be well defined. If QGP and 
hadronic rates are significantly different across T c , ap- 
preciable uncertainties in the calculated photon spectra 
have been found as a result of this ambiguity 14 1. 
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FIG. 1: (Color online) Empirical fits to the direct-photon 
spectrum measured by PHENIX [2|| in p-p collisions at 
y/s = 200 GeV; dash-dotted line: power-law fit [2|§], solid 
and dashed line: x t -scaling ansatz [2ij] with and without K 
factor, respectively. 



B. Non-Thermal Sources 



After subtraction of final-state decays, the main source 
other than "thermal" radiation from the interacting 
medium is associated with "primordial" photons pro- 
duced upon first impact of the nuclei via binary NN col- 
lisions. Primordial photon spectra are usually estimated 
from the direct contribution in pp collisions, where no 
significant rcintcractions arc expected. This is supported 
by the generally good agreement of the measured photon 
spectra with next-to-leading order (NLO) perturbativc 
QCD (pQCD) calculations _26| for primordial produc- 
tion. In the following we will adopt a simple power-law 
fit performed by the PHENIX collaboration @, [27| to 
their pp data 0, S3, cf. Fig. [TJ Fi gure 4 in Ref. [7f re- 
veals that the NLO pQCD calculations are slightly below 
the PHENIX parameterization in the relevant qt range of 
~ 1-7 GeV. Therefore, we will alternatively estimate the 
primordial contribution via a xt-scaling motivated pa- 
rameterization [2^|, upscaled by a if- factor of 2.5 to best 
match the PHENIX pp data for q t = 5-8 GeV. This 
fit produces a slightly smaller yield of photons with mo- 
menta qt = 1-5 GeV compared to the PHENIX fit and is 
thus very similar to the NLO pQCD calculations. 

In principle, photons can also be radiated off fast- 
moving partons (jets) interacting with the medium. The 
combined elli ptic flow of this source is expected to be 
close to zero We therefore subsume this source in 

our primordial contribution (whose fragmentation part, 
e.g., is expected to be reduced in the medium), which 
is a posteriori justified by an adequate description of the 
spectral yields in Au-Au collisions once thermal radiation 
is implemented. 



III. FIREBALL AND TRANSVERSE-FLOW 
FIELD 

The continuous emission of photons throughout the 
evolution of a heavy-ion reaction causes their elliptic-flow 
signal to be more sensitive to its time evolution than that 
of hadronic final states. Therefore, a calculation of the 
photon-«2 requires special care in constructing a realistic 
time evolution of both radial and elliptic flow (the for- 
mer affects the spectral weight of photon emission at a 
given time snapshot). In principle, hydrodynamic models 
are believed to be able to accomplish such a task; how- 
ever, current uncertainties including initial conditions 
(e.g., fluctuations and initial flow fields), viscosity correc- 
tions and the coupling to hadronic cascades in the dilute 
stages render this a challenging task, which has not been 
completed (yet). In the present paper we take a more 
pragmatic (and simple) approach which nevertheless ac- 
curately captures two experimentally established snap- 
shots of the fireball evolution, namely the flow fields at 
chemical freezeout at T ch — 170 MeV and kinetic freeze- 
out at Tf ~ 100 MeV. The latter is well determined 
by the transverse- momentum (pt) sp ectra and elliptic 
flow of light hadrons (n, K, p) |30l - [32l |. while the former 
can be extracted from spectra and flow of multi-strange 
hadrons (lrj l3ll l33l| . A detailed fit of these snapshots us- 
ing an elliptic blast-wave source has been performed in 
Ref. [l5[, which was also shown to be compatible with 
the empirical constituent-quark number and transverse 
kinetic-energy (KEt) scaling of the elliptic flow of light 
and strange hadrons. We use these results to improve 
a previously constructed expanding elliptic fireball 34 1 
so that its evolution passes through these benchmarks 
at the end of the mixed phase and at thermal freezeout. 
Representative examples of the resulting multi-strange (</> 
mesons) and light-hadron (n, p) spectra and t>2 at chemi- 
cal and thermal freezeout, respectively, are illustrated in 

Fig. m 

Concerning initial conditions, we assume an initial lon- 
gitudinal fireball size of zq = 0.6 fm corresponding to 
ct ~ zo/Ay ~ 0.33 fm (Ay ~ 1.8) as our default value 
for both 0-20% and 20-40% centrality classes. With to- 
tal entropies of S = 7900 and 3600 (assumed to be 
conserved), this translates into (average) initial tempe- 
ratures of To ~ 355 MeV and 325 MeV, and charged- 
hadron multiplicities of dN^/dy ~ 610 and 280, respec- 
tively, adjusted to recent STAR data (HJ. 

The time evolution of the inclusive elliptic flow is 
shown in Fig. [3] The observed KEt scaling of light and 
multi-strange hadrons is the main constraint which re- 
quires a rather rapid increase of the bulk v-i, and a sub- 
sequent leveling off shortly after chemical freezeout. The 
final values of ^2.5% and ~5% for the pion-t^ are ad- 



justed to experimental data [31 



Another important aspect in our medium evolution is 
the implementation of chemical freezeout, i.e., the use of 
effective meson (and baryon) chemical potentials to pre- 
serve the observed hadron ratios in the fireball cxpan- 
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FIG. 2: (Color online) Snapshots of p t spectra and t>2 for pions and protons (upper panels), as well as <fi mesons (lower 
panels), following from our fireball evolution in 0-20% Au-Au(y's = 200 AGeV) collisions at thermal and chemical freezeout, 
respectively. The it and p curves are for direct emission only (no resonance feeddown) with absolute normalization while the <j> 
yield is (re-) normalized to the data. Data are from Refs. [lrl l43l. |44| . 
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FIG. 3: (Color online) Time evolution of the inclusive elliptic 
flow for 0-20% and 20-40% Au-Au(^/s=200 AGeV) collisions 
within our fireball model, evaluated with either constituent- 
quark or pion content of the medium. 



sion at temperatures between T c h and Tf . We do this as 
described in Ref. [35j], which was adopted in our previ- 
ous work Q. Most of the hydrodynamic evolutions used 
for photon calculations at RHIC to date assume chem- 
ical equilibrium throughout the hadronic phase. This 
assumption likely leads to an appreciable underestimate 
of the thermal hadronic component in the observed pho- 
ton spectra, and thus of its contribution to the direct- 
photon elliptic flow. For example, typical meson anni- 
hilation processes such as n + p — > 7r + 7 (proceeding 
through t- and s-channel 7r, uj and a\ exchanges), are 
augmented by an initial pion fugacity, z\ = exp(3/z w /T) 
(in Boltzmann approximation), where fx v ~ 100 MeV in 
the vicinity of thermal freezeout, Tf Q ~ 100 MeV. This 
implies a significantly larger enhancement in photon pro- 
duction in the later hadronic stages relative to the con- 
servation of the hadron ratios for which the chemical po- 
tentials are introduced. In other words, the faster cool- 
ing of the fireball in chemical off-equilibrium relative to 
the equilibrium evolution is overcompensated in the lead- 
ing photon-production channels due to a "high" power of 
pion densities. 
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IV. DIRECT-PHOTON SPECTRA 

We start the comparison of our theoretical calculations 
of direct photons to data at RHIC with the absolute 
yields in the transverse- momentum (q t spectra). Let us 
first illustrate the quantitative effect of updating the ra- 
dial expansion starting from our original predictions in 
Ref. [§]. 3 In the latter, a transverse acceleration of the 
fireball surface of ay = 0.053c 2 /fm had been assumed, 
which, together with a fireball lifetime of 15 fm/c, leads 
to a surface velocity of f3 s ~ 0.62 and a freezeout tem- 
perature of Tf = 108 McV for Au-Au collisions in the 
0-20% centrality bin (A part = 280 and iV coU = 765). The 
pertinent photon spectra, displayed in the upper panel of 
Fig- HI closely resemble the results of Fig. 12 in Ref. Q. 4 
A window of QGP-radiation dominance is present for 
q t ~ 1.5-3 GcV. 

The situation changes somewhat with an update per- 
formed in 2007 triggered by the analysis of NA60 dilep- 
tons at the SPS, specifically in the context of their q t - 
spectra. The fireball acceleration was increased to ax = 
0.08-0. lc 2 /fm to better reproduce hadron spectra, which 
also allowed for a significantly improved description of 
the slope parameters in the dilepton q t spectra [l(|. It 
was also checked that the agreement with the WA98 
direct-photon spectra at SPS (36|, as found in Ref. 
was not distorted (see, e.g., Fig. 23 in Ref. [J). At RHIC, 
the pertinent fireball of lifetime r ~ 15 fm/c results in 
a freezeout temperature of Tf Q = 98 McV with a surface 
transverse flow of f3 s = 0.77. The consequences for the 
direct-photon spectra, using the same thermal emission 
rates and fireball chemistry as before, are illustrated in 
the middle panel of Fig. @] while the spectral distribution 
of the QGP radiation is barely affected, the hadronic ra- 
diation spectrum becomes noticeably harder, thus shift- 
ing the crossing with the QGP part up to qt — 1.8 GeV. 
In combination with an improved estimate of the pri- 
mordial emission, adjusted to then available PHENIX pp 
data, the QGP window shrinks appreciably, with a max- 
imum fraction of ca. 42% of the total at qt — 2.1 GeV. 

Finally, recent systematic analyses of light-hadron 



3 For simplicity we will use for this purpose a cylindrically sym- 
metric fireball (no V2) and apply an average boost of 70% of 
the fireball surface flow to the photon spectra in the rest frame, 
</3) = 0.7&. 

4 We note that in Figs. 12 and 13 of Ref. the contribution 
labeled "Hadron Gas" only includes the in-medium p spectral 
function part, not the meson-gas contributions also calculated 
in there. Unfortunately, we recently realized that the spectral 
function part in the photon spectra at RHIC and LHC (Figs. 12 
and 13 in Ref. [ij) was computed with the spin-averaged p prop- 
agator, Dp = (2D T 4- Dp )/3, which, at the photon point (where 
the transverse part, D J , should be used), is by a factor of 2/3 
too small. It was done correctly in the rate plots and for the SPS 
calculations shown in Ref. 0. In the present work we refer to 
the "Hadron Gas" emission as the sum of spectral-function and 
meson-gas contributions. 




q t [GeV] 

FIG. 4: (Color online) The impact of an increasingly strong 
radial flow in an expanding fireball model on direct-photon 
spectra in 0-20% central Au-Au collisions at RHIC. The 
transverse fireball acceleration increases from 0.053/fm (up- 
per panel, corresponding to Ref. [jj) via 0.08 (middle panel; 
see Refs. [3, G]}) to 0.12/fm (lower panel). The same QGP 
and hadronic emission rates have been used in all cases, while 
the primordial contribution has been upscaled in the middle 
and lower panel. As a benchmark, we also show the pertinent 
PHENIX data 0. 



spectra by the STAR collaboration [32j requires an even 
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FIG. 5: (Color online) Comparison of our calculated direct-photon spectra (left panels) and their elliptic-flow coefficient (right 
panels) from an elhptically expanding fireball model with QGP and hadronic radiation, supplemented with primordial emission, 
to PHENIX data @,[H in 0-20% (upper panels) and 20-40% (lower panels) central Au-Au(Vs = 200 AGeV) collisions. Models 
(a) and (b) in the right panels refer to the use of the pQCD parameterization and the PHENIX fit for primordial production, 
respectively (in the left panels, only model (a) is displayed). 



harder expansion, to reach a thermal freezeout configu- 
ration at Xfo ~ 90-95 McV and (/3) ~ 0.59 for central 
Au-Au collisions. This can be achieved in the fireball 
model by a further increase of the acceleration to ot = 
0.12 c 2 /fm at a slightly reduced lifetime of r ~ 14 fm/c. 
As expected, for the direct-photon spectra this implies 
a further hardening of the hadronic emission and thus 
an additional squeezing of the QGP window to a small 
region around qt ~ 2.4 GeV, cf. lower panel of Fig. |4] 
The increasing transverse flow of the three fireballs also 
seems to improve the description of the PHENIX direct- 
photon spectra, although that was not the objective of 
this exercise. 

The final (third) fireball setup has been refined by im- 
plementing realistic ellipticities and an explicit linearly 
increasing flow profile in the transverse boost of the pho- 
ton emission rate. The corresponding comparisons to 
hadronic data have been discussed in the previous sec- 
tion. The resulting direct-photon q t spectra from thermal 



QGP and hadronic sources, supplemented with an ^en- 
sealed primordial contribution, are compared to PHENIX 
data in the 0-20% and 20-40% centrality classes of Au- 
Au(y/s = 200AGeV) in the two left panels of Fig. [5] 
The more central data set is fairly well reproduced, even 
though there appears to be a slight underestimation of 
the datum at the lowest q t ~ 1.2 GeV. The inclusion 
of the full transverse-flow profile leads to a further hard- 
ening of the hadronic component relative to the lower 
panel in Fig. [U while the QGP component is essentially 
unaffected, even at the highest qt (continuing the con- 
stant trend of the three panels in Fig. [4j. This means 
that the high-qt QGP radiation is entirely determined by 
the earliest radiation, where no flow has built up yet; the 
subsequent QGP flow cannot overcome the softening due 
to the decreasing temperature. This further implies a 
significant dependence of the high-q t QGP yield on the 
thcrmalization time (a quantitative example will be dis- 
cussed at the end of Sec.fV]). On the other hand, hadronic 
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p T (GeV) p T (GeV) 

FIG. 6: (Color online) Same as Fig. [S] for 0-20% Au-Au collisions, but with a QGP contribution evaluated for a reduced 
thermalization time of tq ~ 0.17 fm/c translating into an average initial temperature of To ~ 445 MeV. 



emission only sets in at T c when there is already substan- 
tial flow in the system, and thus even at high momenta 
the hadronic spectra arc sensitive to the fireball flow field. 

In the 20-40% centrality bin, the discrepancy between 
the theoretical yields and the data becomes somewhat 
more severe, hinting at a missing relatively soft source 
(and therefore suggestive for the later hadronic phase). 
One speculation at this point could be related to u — > n°j 
decays. These have been subtracted by the PHENIX col- 
laboration employing m t scaling of the uj spectra with 



27. |, assuming lo/tt = 1, as found in pp mea- 
surements [22l, as well as in 0-92% Au-Au collisions for 
pt > 4 GeV. //, however, u> mesons at lower pt become 
part of the chemically equilibrated medium in heavy-ion 
collisions, one expects their multiplicity at given m t to 
be up to 3 times larger, due their spin degeneracy. In 
this case there might be a direct-photon component in 
the Au-Au data at low q t < 2 GeV due to some frac- 
tion of final-state u> — > tt°j decays which have not been 
subtracted (and which would carry large v%). This possi- 
bility may be worth further experimental and theoretical 
study. 

It is quite remarkable that the hadronic yield domi- 
nates over the QGP one over the entire plotted range. 
This will have obvious ramifications for the 1*2 of the 
direct photons, which is larger in the hadronic phase. 
The sub-leading role of the (early) QGP component fur- 
ther implies that the effects of initial-state fluctuations 
on thermal-photon production [20L l2ll ] are diminished. 

To examine the dependence of the QGP yield on the 
thermalization time, we have conducted calculations with 
a factor-2 reduced initial longitudinal size, zq = 0.3 fm, 
corresponding to tq ~ 0.17 fm/c as used, e.g., in Ref. [l4j |. 
cf. Fig. © The QGP spectra in 0-20% Au-Au collisions 
increase over the zq = 0.6 fm calculation by a factor of 
1.6, 2.7 and 4.8 at q t = 2, 3 and 4 GeV, respectively, and 
turn out to be in fair agreement (within ca. 30%) with 
the hydrodynamic calculations reported in Ref. [14[ (us- 




ing smooth initial conditions). The significance of this 
increase mostly pertains to momenta, q t > 2 GeV, where 
a small "QGP window" reopens, but it does not signifi- 
cantly affect the description of the experimental yields. 

To further characterize the nature of the direct-photon 
excess (i.e., beyond the pp-scaled primoridial emission), 
we evaluate the effective slope parameters, T e g, of our 
thermal spectra. We recall that PHENIX extracted the 
effective slope of the excess radiation in their data as 
T cff = 221±19 stat ±19 syst MeV Q. In Fig. [7] we compare 
this range with the temperature evolution, T(t), of our 
fireball; they only overlap inside the QGP phase. How- 
ever, when accounting for the flow-induced blue shift, as 
estimated by the schematic expression for a massless par- 
ticle. 



(1) 



the overlap with the experimental window is shifted to 
significantly later in the evolution, mostly for a flow- 
ing hadronic source with a restframe temperature of 
T ~ 100-150 MeV. This suggests a rcintcrpretation 
of the experimental slope as mainly hadronic in origin, 
which, as we will see in Sec. |V] below, is further sup- 
ported by the V2 data. An explicit fit of the slope to our 
total thermal spectrum from the elliptic fireball (with 
T = 355 MeV) in the range q t ~ 1 - 3 GeV yields 
T c ff ~ 240-250 MeV, which is at the upper end of the data 
(consistent with the slight underestimate of the lowest-gt 
datum; also note that the use of the average, ((3) = 0.7/3 s 
in Eq. ([1]), tends to underestimate the actual slopes, es- 
pecially at high qt and j3 s ; we noted that already when 
going from the spectra in the lower panel of Fig. @] to 
the full results in the upper left panel of Fig. [5]). Higher 
initial temperatures are less favorable, since they result 
in a further increase of the slope, e.g., by 10-15 MeV for 
T = 445 MeV. 
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FIG. 7: (Color online) Time dependence of the effective slope 
parameter of thermal-photon radiation emitted from our fire- 
ball for 0-20% Au-Au collisions (solid (blue) line, given by 
QGP and hadronic sources for r — to < 3 fm/c and > 6fm/c, 
respectively), using Eq. ([TJ with (/?) = 0.7/3 s . For comparison, 
we plot the effective slope extracted by PHENIX from their 
data (horizontal line with dashed lines indicating the experi- 
mental error) Q, and the "true" temperature in the thermal 
restframe (long-dashed line). Also shown is the surface flow 
velocity of the fireball (dash-dotted line) . 



V. DIRECT-PHOTON ELLIPTIC FLOW 

With a fair description of the photon qt spectra at 
hand, we proceed to the calculation of the direct-photon 
elliptic flow as a function of transverse momentum, 
v 2 (<?*)• The results for the 0-20% and 20-40% centrality 
classes of Au-Au collisions at RHIC arc again compared 
to PHENIX data, see upper- and lower-right panel of 
Fig. [51 respectively. The shape of our calculated ^2(9*) 
matches the measurements rather well, but our maxi- 
mum value of ca. 7.5% (12%) is below the central values 
of the 0-20% (20-40%) data. Our results are a factor of 
3-4 larger than existing calculations using hydrodynamic 
expansions, and reach into the lower end of the (mostly 
systematic) experimental uncertainties. The main differ- 
ences compared to the hydro calculations are the follow- 
ing: our equilibrium hadronic rates Q are significantly 
larger than in previous studies, our hadronic phase in- 
cludes meson-chemical potentials and lasts longer (both 
due to a smaller Tf as dictated by data and a slightly 
larger T c h as used in our previous calculations [1, Qjl, EDI), 
and our bulk elliptic flow is built up faster than in 
standard (ideal) hydro calculations (recall that only a 
small increase of the V2 during the hadronic phase fa- 
cilitates its KEr-scaling properties of multi-strange and 
light hadrons [15j ; also note that a more rapid expansion 
in the QGP and transition region has been identified as 
an important ingredient to solve the "HBT puzzle" EDI ) ■ 
The combined effect of these four points is a thermal 
source which is dominated by hadronic emission carrying 
most of the finally observed elliptic flow from its begin- 
ning on, i.e., for T < T c . For the reduced thermalization 



time of To = 0.17 fm/c the maximal vi drops by ca. 15% 
to v? 3 * ~ 6.3%, cf. right panel in Fig.© 

In the present work, we have assumed a critical temper- 
ature of T c ~ 180 MeV, which is in line with N f = 2 + 1 

flavor lattice calculations reported, e.g., in Ref. [4(ij |. 

However, very recent lattice data [4l| (as well as ear- 
lier ones [42|) indicate that T c could be as low as 155- 
160 MeV and therefore one should ask what impact this 
could have on our results. A pertinent study has been 
done in the context of dilcpton production at SPS en- 
ergies (Toj . When varying the critical temperature by 
±15 MeV around the default value of T c = 175 MeV, the 
QGP emission spectra in the intermediate-mass region 
(M > 1 GeV) vary by up to ±50% (less at higher masses 
where the contribution from earlier phases increases) . At 
the same time, the hadronic emission part varies by ap- 
proximately the same amount in the opposite direction, 
so that the total yield roughly stays the same while the 
relative QGP and hadronic partition varies appreciably. 
We expect similar effects for the photon q t spectra. If one 
still requires the multi-strange hadron spectra to freeze 
out at T c , we expect that the V2 in the hadronic phase 
does not significantly change with T c . However, since for 
smaller T c the hadronic contribution to the direct pho- 
ton yields is reduced, so should be the total (weighted) 
V2- For example, if for the 0-20% centrality class the 
QGP— hadronic partition of 1/3—2/3 at q t ~ 2 GeV 
changes to 1/2—1/2, we estimate that i>2(<?t) is lowered 
from 7% to 6%. 



VI. CONCLUSIONS 

We have updated and extended our calculations of 
thermal-photon spectra at RHIC by constructing an im- 
proved elliptic fireball expansion which is quantitatively 
constrained by bulk-hadron data. In particular, we have 
implemented the notion of sequential freezeout by repro- 
ducing an empirical extraction of radial and elliptic flow 
from multi-strange and light-hadron data at chemical and 
kinetic decoupling, respectively. The fireball evolution 
has been combined with existing photon emission rates in 
the hadronic and QGP phases to obtain thermal-photon 
spectra in Au-Au collisions at RHIC. Supplemented with 
a primordial component estimated from pp collisions, we 
have compared our calculations to spectra and v 2 of di- 
rect photons as recently measured by the PHENIX col- 
laboration. Due to a large medium flow (as required 
by hadron data), relatively large hadronic photon rates 
(approximately degenerate with QGP rates around T c ), 
and effective chemical potentials to conserve the observed 
hadron ratios, we have found that the hadronic medium 
outshines the QGP for most of the momenta where ther- 
mal radiation is relevant. This, in turn, leads to a ma- 
ximal elliptic flow coefficient of V2 — 10% in semicentral 
Au-Au, which is a factor of ~ 3 increased over previ- 
ous estimates based on QGP-dominated emission. Con- 
sequently, the discrepancy with the PHENIX V2 data is 
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reduced appreciably. Our results arc corroborated by 
evaluating the effective slope parameters of the radiation 
from the thermal source, which have the largest overlap 
with the experimental value of T ff — 220 MeV in the 
(flowing) hadronic phase. Initial QGP temperatures of 
above To ~ 400 MeV are increasingly disfavored by both 
the slope and vi data. Further scrutiny is needed whether 
these results can be confirmed in dynamical space-time 
models, i.e., in hydrodynamical and transport simula- 
tions, and what the quantitative impact of the transition 
temperature is (u sing , e.g., the (3+l)D viscous hydro- 
dynamics of Ref. [45|). The prevalence of the hadronic 
emission in our calculations reiterates the necessity of 
a good understanding of the strongly coupled hadronic 
phase in heavy-ion collisions. With these considerations, 
a satisfactory explanation of the (surprisingly?) strong 
direct-photon vi signal at RHIC might be possible. The 



extension of the studies to virtual photons (aka dilcp- 
tons) would also be illuminating. First phenomenological 
studies of this observable have been initiated [46|, H3] and 
are also planned within our framework. 
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